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Abstract 

Neutron-deuteron scattering in the context of "pion-less" Effective Field Theory 
at very low energies is investigated to next-to-next-to- leading order. Convergence is 
improved by fitting the two-nucleon contact interactions to the tail of the deuteron 
wave-function, a procedure known as Z-parameterisation and extended here to the 
three-nucleon system. The improvement is particularly striking in the doublet-S wave 
(triton) channel, where better agreement to potential-model calculations and better 
convergence from order to order in the power counting is achieved for momenta as 
high as ~ 120 MeV. Investigating the cut-off dependence of the phase-shifts, one 
confirms numerically the analytical finding that the first momentum-dependent three- 
body force enters at N^LO. The other partial waves converge also substantially faster. 
Effective-range parameters of the nd-system are determined, e.g. for the quartet-S- 
wave scattering length aq = [6.35 it 0.02] fm, which compares favourably both in 
magnitude and uncertainty with recent high-precision potential-model determinations. 
Differential cross-sections up to E'lab ~ 15 MeV agree with data. 
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1 Introduction 



The few-nucleon system at very low energies is an important tool to understand key physical 
questions. On the one hand, a new generation of high-precision experiments at the lower 
end of the energy spectrum use both polarised targets and beams, and employ neutrino 
detectors or radioactive-beam facilities to extract e.g. neutron properties, neutrino masses 
and reactions relevant for nuclear astro-physics. On the other hand, a plethora of pivotal 
physical processes is hard to access directly in experiments, like reaction rates in Big-Bang 
nucleo-synthesis. In both cases, it is mandatory that binding effects in light nuclei are taken 
into account with as few bias as possible towards a particular model of the few-nucleon 
system. Theory must provide such model-independent extraction and calibration methods 
or predictions at very low energies, often well below 20 MeV. 

While precision - the numerical stability of computations - can be controlled by combin- 
ing sophisticated algorithms with Moore's law and is hence (albeit sometimes formidable) 
not a fundamental problem, the necessary accuracy - namely an estimate of the system- 
atic uncertainties of a theoretical ansatz - is harder to win and mandates understanding the 
physical system at hand. In the three-nucleon sector, this calls in particular for a systematic 
understanding of the role three-nucleon forces play. Traditionally, these were often intro- 
duced a posteriori to cure discrepancies between data and calculations, e.g. for the triton 
binding energy, but such a path is of course untenable when data are scarce or absent. There 
is also a host of deviations between experiment and theory for which a suitable three-body 
force could not yet be constructed, like the famed A^-problem [2]. 

The so-called "pion-less" version of Effective Field Theory in Nuclear Physics (EFT(^)) 
aspires to provide just such a systematic classification of all forces. At its heart lies the 
tenet that Physics at those very low energies can be described by point-like interactions be- 
tween nucleons only: One cannot identify pions as the lightest exchange-particles between 
nucleons as long as the typical external momentum ptyp in a reaction is below the pion mass 

because the Compton wave-lengths are not small enough to resolve the nuclear forces 
as originating in part from one-pion exchange. With all particles but the nucleons thus 
"integrated out" as heavy, one can identify a small, dimension-less parameter Q = -C 1, 

where ~ is the typical momentum scale at which the one-pion exchange is resolved 
and EFT{f) must break down. The resultant power-counting orders each process according 
to the power in Q at which it starts to contribute and establishes therefore an ordering 
scheme which is used for systematically improvable, rigorous error estimates. This and 
the systematic, gauge invariant inclusion of external electro-weak currents and relativis- 
tic effects, see e.g. distinguishes EFT(^) from its historical roots, the Effective Range 
Expansion and the model-independent approach to three-body physics 

EFT(^) is a mathematically well-defined, systematic low-energy theory of QCD. It is 
computationally considerably simpler than potential models or the "pion-ful" version of 
EFT, which attempts to extend Chiral Perturbation Theory to the few-nucleon system [HI 13 
IH]. Conceptually, many problems which are also found when formulating a fully consistent 
"pion-ful" EFT are encountered in a simpler setting. It was used quite successfully to 
provide model-independent results for a cornucopia of two-nucleon processes with external 
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electro-weak probes, see the reviews El E] and references therein, also for a sketch 
of its development. Bedaque et al. ^21 showed that a momentum-independent three- 
body force must be LO for the triton and sketched the path to include further three-body 
forces systematically at higher orders. Recently, a rigorous power-counting to all orders 
was developed for the three-nucleon forces in this approach, opening the path to new high- 
accuracy extractions and predictions of nucleonic and nuclear properties ^1]. Barford and 
Birse confirmed this in an analysis of the renormalisation group flow of the position-space 
version of the problem ^Hl- Now, convergence issues become interesting. They are the 
subject of this article. 

In EFT(^), the strengthes of the contact interactions between the nucleons can be deter- 
mined from low-energy observables in various ways, which differ in principle only by higher- 
order effects. Still, it is standard practise to to improve the speed of convergence by physical 
considerations. One would for example not try to describe deuteron properties by starting 
from the Effective Range Expansion (ERE) in the '^Si-channel of iVA^-scattering around 
zero momentum, in which the correct position of the deuteron pole at Bt = 2.2246 . . . MeV 
is reached only perturbatively. Instead, it is prudent to put the deuteron pole-position in 
the right place, and give the pole its correct strength. The deuteron wave-function has then 
the proper asymptotic fall-off and normalisation. Phillips et al. ^Hl showed that this Z- 
parameterisation is an effective way to sum up the dominant effective-range contributions. 
It both simplifies calculations of deuteron properties and improves convergence. 

As the ^Si-channel of A^A^-scattering is an important sub-cluster of the three-nucleon 
system, it is natural to consider the implications of this choice for the triton and other 
partial waves. With a slight extension, it can be included into the standard EFT treatment 
of the three-nucleon system in which an auxiliary two-nucleon field is introduced to simplify 
computations. This article will show that not only is the convergence both from order to 
order in the expansion and to experimental phase-shifts greatly improved; its results support 
also the power counting for the three-body forces pertinent to the triton - mathematically 
rigorously proven to all orders in ^3] - by an error analysis of the ^Si-channel. 

This article is organised as follows: The next Section merges Z-parameterisation with the 
auxiliary-field method. After a brief review of the formalism to compute Nd scattering 
in EFT(^) and the occurrence of three-body forces in the triton channel at the beginning 
of Section El the results for the ^Si-channel, triton properties and the Phillips line are 

_____ 2 

examined in Sub-Sect. 13.31 I touch upon the pros and cons of including the effective range 
as leading-order in Sect. 13.41 Finally, the effect of Z-parameterisation on the other partial 
waves of the A^rf-system is discussed in Sub-Sect. 13.51 Various effective-range parameters 
of low partial waves are also listed, together with comparisons between experimental and 
theoretical cross-sections. The Conclusions in Sect. I^are followed by an Appendix on the 
derivation of the three-body equations in EFT(|f), defining also the spin-isospin projection 
operators in the three-nucleon system. 
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2 Z-Parameterisation and Auxiliary-Field-Formalism 
in the Two-Nucleon System 

2.1 Merging the Auxiliary Field Formalism . . . 

It is a standard technique to simplify calculations in EFT(^) both in the two-and three- 
nucleon system by considering not directly contact interactions between four nucleons, but 
to introduce auxiliary fields with the quantum numbers of the two-nucleon real and virtual 
bound states, coupling to two nucleons ^, see e.g. [T7| ITHl IT^ QUI 1^: 

Atv,* = -y\(fhN^PiN)+YL.c] (2.1) 

n+l" 




pp. / f)2 ^2 

A*-coJ 1^0 + ^ + 5- -Vc„Ji9o + -^+^* 



AM M ^ \ " AM M 



Here, = (^) is the nucleon iso-doublet and = -^T20-2cr^ the projector onto the spin- 
triplet iso-spin-singlet state with vector index 2 = 1,2, 3. a* (r^^) are the spin (iso-spin) Pauli 
matrices, A = 1,2,3 the iso- vector indices. The auxiliary field dt represents the deuteron 
with binding momentum 7^ = y/MB^ = 45.7025 MeV. Strictly speaking, one should replace 
the (iso-scalar) nucleon mass M — M — Sf/2 to obtain the correct deuteron mass, but this 
effect is negligible in what follows. Only the ^Si-channel of NN scattering is considered 
here; the ^So-channel is discussed in Sect. 12.31 The choice of sign for Cnt is traditional [T7] . 

That this Lagrangean is on-shell equivalent to the one containing only nucleon fields 
was formally shown in [2111 by a GauB'ian integration over dt, followed by a field-redefinition 
and disregarding terms with more than four nucleon fields. The advantage of this scheme 
- stressed repeatedly, e.g. in (TU [THJ UHl I2D1 ^ is that the parameters follow naive 
dimensional analysis: As again discussed in Sect. 12.21 ~ Q, with mass-dimension 1, is 
LO; the dimension-less parameter cot ~ first appears at NLO since it comes with two 
powers of momentum, CotP^ ~ Q^', while a dimension-ful operator proportional to Cnt ~ 
enters at N^"'"'"^LO because it is accompanied by 2n + 2 powers of the typical momentum. 

All interactions which are not pure S-wave - like SD-mixing, P-wave scattering between 
two nucleons etc. - are added either as interactions between the deuteron dt and two nucle- 
ons, or between four nucleons [211 122] • They are not listed as they are of higher order than 
necessary in the following - with the exception of SD mixing, whose effect will be neglected 
in the three-nucleon system, see the discussion in Sect. 13.11 The only additional term at 
N^LO is the kinetic energy of the nucleon, which up to relativistic corrections is: 

£1^ = N^{ido + ^]N (2.2) 



^As by-product, I attempt to unify the notational cornucopia of which the present author is not com- 
pletely innocent. Throughout, the sub-script t (s) denotes quantities in the spin-triplet (singlet) channel of 
NN scattering. 
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The bare deuteron propagator i/A^ is thus dressed at LO by all interactions proportional 
to y with an arbitrary number of loops, see Fig. ^ At NLO, one perturbative insertion 
proportional to Cgt is included, followed by two at N^LO, and in general n at N"LO. In 
addition, one insertion of the operator proportional to Cnt enters at N^^+^LO, etc. 



+ o + 



)( )( )( 



icit(po-fe-+ 



Figure 1: Top: Re- summation of the bare deuteron propagator (thick gray line) into the 
dressed deuteron propagator (double line) at LO by dressing with two-nucleon bubbles. 
Bottom: The NLO, N^LO and N^LO corrections to the deuteron propagator. The cross 
denotes an insertion of the deuteron kinetic-energy operator proportional to Cot, the star 
one proportional to Cu- 



Choosing for the dtNN coupling constant without loss of generality 



47r 

M 



(2.3) 



the auxiliary- field propagator with kinetic energy pq and momentum p becomes particularly 
simple in the low-energy expansion: 



i Dt{po,p) 



'^t + - \/ T - ^Po - 1^ "^=0 



n=l ^ 



n+l 



:2.4) 



- Mpo - ie - cotipo - |g + 17) - E Cnt - |g + 

n=l ^ 



M 



Here, /x is the regulator of the linear divergence in the nucleon loop, regularised using 
dimensional regularisation with the PDS subtraction scheme [22] • The scattering amplitude 
between two non-relativistic nucleons with relative centre-of-mass (cm) momentum k in the 
^Si-channel is obtained by multiplying with —y"^ = —An/M and setting the nucleons on-shell 
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M At + fx + 



An 



ik ^ 

m=0 



I i.2 



M' 



CO 

n=l 



fci _L I? 

M Af 



At + /i + i/c 



Aj + /i + ifc 



oo / 

E c„t ( 

n=l ^ 



A/ Af 



n+1' 



:2.5) 



To sum all effective-range corrections to all orders in the second line of ()2.4|1 and (j2.5|l has no 
advantage but to shorten the following determination of the low-energy coefficients; actual 
calculations involve only a finite number of c„t's, in which case the difference between the two 
lines is formally of higher order. Problems with re-summing the effective-range corrections 
"to all orders" are discussed in Sect. 13.41 



2.2 ... with Z-Parameterisation 



In Effective Range Expansion (ERE) around the deuteron pole jl] , the scattering amplitude 
of two non-relativistic nucleons with relative cm momentum k in the ^Si-channel is 



^7V7v(^Si,/i;) 



M 
Air 



[2.6) 



k^)-T.Pnt{lI + k^Y+^ + ik 

n=l 



where pot = 1-764 fm is the effective range, pu = 0.389 fm'^ the shape parameter etc. 
Matching ()2.5|) to the ERE-result ()2.6|) . one obtains the ERE-parameterisation of EFT(^): 



At,ERE + A'ERE — 7i ! Cot,EKE 



Pot 



M , c„t,ERE = Put M"+i Vn > 1 



(2.7) 



While 7t is clearly a typical low- momentum scale in EFT(|f), 7^ ^ ~ m^, the other 
parameters encode Physics beyond the breakdown scale and hence one expects pnt ~ A^^"^^ 
from dimensional analysis. Indeed, contributions to the scattering amplitude for example 
from the shape-parameter are small, with the expansion parameter Q estimated as Q ~ 
ItPit^ ~ 0.17. Therefore, the amplitude can be expanded in powers of Q ~ ll^^'^'^'^^^ Pnt, and 
an error estimate follows, rendering Effective Range Theory useful: 



ANN^Suk) = - 



M 



1 



An 'jt + ik 



Pot It +k'^ , f Pot 7? + k 



2 -ft + ik \ 2 -ft + ik 



(2.8) 



A problem with the expanded version ()2.8|) arises however because the effective range is 
numerically somewhat larger than naively expected, po > "^^^ or 



-ftpot ~ 0.41 



(2.9) 
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While the deuteron binding energy is in ERE-parameterisation reproduced immediately at 
LO, the residue of the scattering amplitude at the deuteron pole is now only expanded in 
this numerically not-so-small parameter, 

Res[^ArAr(^Si, k = I'jt)] = -j^ Zf with 

= 1.690(3) (2.10) 



1 - ItPot 

= I + Itpot + iltpot)^ + iltpotf + . . . = 1 + 0.409 + 0.167 + 0.068 + ... , 

so that even including N^LO-effects, seven percent are missed to close the gap to the exact 
value. In contrast, not considering the N^LO shape-parameter effects in A leads only to a 
deviation on the order of "yfpu = 0.5%. 

Clearly, the residue is an important characteristic of the two-body system as it determines 
the asymptotic normalisation of the deuteron wave-function at large distances r: 



^'dcutcron(r ^ Oo) = W . (2.11) 

Phillips et al. ^Hj re-summed therefore the series in 'jtPot partially by expressing pot via Zt. 
To discuss its advantage, re-write (j2.6j) in a superficially more complicated way as 



l + {Zt 



^'^"^ 1 + ^ (l + S) - Ztt^Pnt ilt - ik) (7| + F)" 



(2.12) 



^ [1 + (Z, - 1)] 



7t + 



Formally, Zt — 1 ^ 0.69 is still treated as a small expansion parameter in the resulting 
Z-parameterisation, but the deuteron residue is now restored already at NLO with no cor- 
rections from higher orders as the second denominator in ()2.12j) has no residue at '-ft = —ik: 



Zt= 1 +{Zt-l)+ , n>2. (2.13) 

^To^ ^ NLO ' N^^io 

Outside the deuteron pole, the additional terms from Zt also converge faster because the 
second term in the second denominator is a power series in {Zt — l)/2 ^ 0.3. A slight 
dis-advantage is that the expansion parameter for the higher-order correction is now bigger 
by 70%: Zt Pntl't^^^ ■ Thus, the N^LO-correction from the shape parameter is now ~ 0.8%. 
At this level of accuracy, however, other corrections (e.g. from P-wave interactions) must be 
considered, whose scale is also set by the pion mass and which are estimated to be stronger 
because (7f/(A^ ~ ""^Tr))^ ~ 3%. 

In Z-parameterisation, the expanded version of the EFT(^)-amplitude ()2.5p is therefore 
first matched to reproduce the correct deuteron pole position: At + p = 7*. Then, the 
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residue is found starting at NLO as 



n=l 



^. = l + E(^o.^) . (2.14) 



M 



At LO, no additional free parameter exists in EFT and the residue is one. At higher orders, 
the residue is set equal to Zf when cqi is suitably chosen: With the expansion truncated 

at some finite order, Cot = Yl % contains now contributions from higher orders in the 

n=0 

Q-expansion, Cg^^ being 0{Q"). To summarise, the parameters in Z-parameterisation are 

= ^-Q' (2.15) 

At + fi = it--Q (2.16) 

M 

2^ 

Cnt = Put M"+i ~ go Vn > 1 . (2.18) 

In contradistinction, cqi receives in ERE-parameterisation only a contribution at order Q^, 
and none at higher orders ()2.7p . However, such an expansion is with our choice of Lagrangean 
()2.H1 only encountered for one parameter, and not - as previously - for all [T ^ IT? H OI H I22j. 
The auxiliary-field propagator to N^LO, 

Dt{po,P) 



£^ = (-)« (z, - 1)"+^ ~ g'^ (2.17) 



7t- V^-Mpo-ie 



1 + 



LO 



27* 



It 



27* 



£__Mpo-7t )+o(g 



NLO 



N^LO 



(2.19) 



is thus a faster-converging alternative to the usual ERE of the A^A^-scattering amplitude 
in the ^Si-channel. The wave-function renormalisation, i.e. the residue of the deuteron 
propagator, is by construction exact at NLO: 



d 



dpo A(po,0) 



Po=-Tr. 



27* 
M 



1 +{Zt-l)+ , n > 2 

LO NLO N"LO 



(2.20) 



One could also replace directly the deuteron propagator in ERE by the result in Z-parameter- 
isation ()2.12j) . as e.g. argued to be computationally simpler by Beane and Savage [23 • This 
would be mandatory if 7tPo* ~ 1- Equivalently, if one re-sums all orders in Cgt, one finds 
again the ERE-parameterisation ()2.7p . cqj = ~ ^ However, there are a number 

of dis-advantages of this approach, as will be discussed in Sect. 13.41 
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2.3 Z-Parameterisation for the ^So-channel? 

The obvious question is: Why not also impose Z-parameterisation for the ^So-channel of 
NN scattering? As its bound state is only virtual, one usually performs the ERE around 
zero momentum, 



^AfAr(^So, k) — — 



M 
in 



1 



— -"-tk^-H Tnsk^''^^ + i k 



(2.21) 



with as = 23.714 fm, ros = 2.73 fm, ris = —0.48 fm^ the scattering length, effective range 
and shape parameter. As ros/os ~ 0.11, the series-expansion of the residue of the virtual 
bound-state converges much faster than in the deuteron channel. Still, in order to simplify 
notation, consider a Lagrangean analogous to the ^Si-channel: 



J0.2N,s = -y {df {N^P^N) + E.c.) 



(2.22) 



As - Cos ido 



M ^ 



4M M 



n=l 



92 2 

^^°+4M + Af 



The auxiliary field ds represents the spin-singlet iso-spin-triplet state, whose projector is 
= ■^a2T2T^. One now first re- writes the ERE expansion ()2.2H1 in the form analogous 
to ()2.6p in which the pole position does not change from order to order. 



^ArAr(^So, k) 



M 

An 



(2.23) 



7s - ^ (7.2 + k')-E PnM + k^Y^^ + i k 

n=l 

and determines the coefficients by matching as 

00 

1 , ^Os 2 / 1 

- + ^7s - 2^'^n.(-7. 

Ujs ^ -1 

n=l 

00 

ro, + 2 J](n + 1) r^s{-ll 
'n + 1" 



Is 



POs 



,2\n+l 



r'sT 



n=l 

00 



Pv 



''ns[-ls 



2\n-l 



m > 



(2.24) 
(2.25) 
(2.26) 



Truncation at pos (Pis) leads to the numerical values 7^ = —7.8904 MeV (—7.8902 MeV), 
Pos = 2.730 fm (2.733 fm) (, pis = ris). The residue is indeed very close to unity: 

1 



1 - IsPOs 



0.9016 (0.9015) 



(2.27) 



At N^LO, the difference between the perturbatively built residue 1 + 'jsPos + {isPosY = 
1 — 0.1092 + 0.0119 = 0.9027 and the exact values is with 0.1% considerably smaller than 
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leaving out relativistic and other effects. In the results presented in the following Sec- 
tion, the difference between the ERE-parameterisation and Z-parameterisation for the ^Sq- 
channel cannot be discerned in the plots of the phase-shifts. Still, the re-formulation serves 
the purpose to compactify and simplify formulae. The parameters of the Lagrangean are 
determined in Z-parameterisation by the analogue to ()2.15|1 : 

(2.28) 
(2.29) 

— r^Q"" (2.30) 
^ Q° Vn > 1 . (2.31) 

Note that in contradistinction to previous work, y is chosen to be identical in the spin-triplet 
and spin-singlet channel and thus carries no sub-script. 



A, + /i 


M ^ 

= is^ Q 


An) 


= (-) 






Pns 





3 The Three-Body System in Z-Parameterisation 
3.1 Formalism 

With the parameters of the two-nucleon Lagrangean fixed by Z-parameterisation, it is now 
straight-forward to state the equations governing neutron-deuteron scattering. As they were 
derived repeatedly in the literature (see e.g. [131 122 Ell), following presentation focuses 
mainly on notation. Appendix ^ contains a brief overview, defining also the pertinent 
projection operators on the various partial waves and three-nucleon configurations. 

Two cluster-configurations exist in the three-nucleon system: The A^dt-cluster with total 
spin S = I or S* = |, depending on whether the deuteron and nucleon spins are parallel 
or anti-parallel; and the iVcis-cluster which has total spin S — 2 5 

is a scalar. The 

leading-order three-particle amplitude is 0((5~^) (before wave- function renormalisation) and 
includes all diagrams built out of the leading two-body interactions, i.e. the ones proportional 
to y and A^, A^ in the two-nucleon Lagrangeans ()2.1|/I2.2^ . The resultant Faddeev integral 
equation - first derived by Skorniakov and Ter-Martirosian [21] without three-body force - 
is pictorially represented in Fig. |21 

As the Lagrangean up to N^LO does not mix partial waves or flip the spin of the auxiliary 
fields, angular momentum is conserved in the quartet and doublet channels. Strictly speak- 
ing, SD-mixing in the deuteron channel produces a splitting and mixing of the three-body 
amplitudes with the same spin and angular momentum but different total angular momen- 
tum. However, we limit ourselves to the averaged phase-shifts, as the spitting is at this order 
not going to be realistic enough to describe spin-observables in the three-nucleon system. 
This path was also pursued in and [22] which used ERE-parameterisation, to which the 
findings in Z-parameterisation will be compared. An analysis of these spitting is postponed 
to a future presentation, which also deals with spin-observables in the three-nucleon system. 
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Figure 2: The Faddeev equation for A^d-scattering to N^LO. Thick sohd hue: propagator of 
the two intermediate auxihary fields dg and dt, denoted by V, see ()3.6|) : /C: propagator of 
the exchanged nucleon, see ()3.2|) : H: three-body force, see ()3.11|) . 



The spin-quartet channel of the A^d-system is particularly simple. The Zth partial wave of 
the amplitude tq^ in the centre-of-mass (cm) frame is given by ^ 

oo 

ti;\E;k,p) = -A7rJC^'\E;k,p) + ^ J dq q' }C^'\E;q,p) Dt{E - £^,q) t^^\E-k,q) .(3.1) 



The total non-relativistic energy is E := fij — the incoming (outgoing) deuteron- 

momentum k (p); and the projected propagator of the exchanged nucleon on angular mo- 
mentum / ^ 



}C^^\E;q,p) := - / dr 



Pi(x) 



p2 _|_ g2 _ ]\/[jt; — ie + pqx 



{-ly ^ [p^ + g2 _ ME - ie 
Qi 



pq 



pq 



(3.2) 

where the /th Legendre polynomial of the second kind with complex argument is as in [221 



Qi{z) 



m 

z-t 



(3.3) 



In the doublet channel, the Faddeev equation is two-dimensional in cluster-configuration 
space as both Ndt- and A^tig-configurations contribute: 



^\E;k,p) = 27r 



/C«(i?;fc,p) 



1 



5'° n{E; A) 



1 



(3.4) 



A) 



1 -1 

-1 1 



^V{E-^,q)^^{E-k,q) 



^The sub-script q (d) denotes quantities in the spin-quartet (doublet) channel of the three-nucleon system. 
^This corrects a typographical error in a previous paper on the higher partial waves in the A'^d-system |22| . 



affecting the odd partial waves. 
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The three-body force Ti. will be discussed in the next Sub-section. The vector 



4') := ( 5f ) (3.5) 

is built out of the two amplitudes which get mixed: for the Ndt Ndt-piocess, and 
t^j\^ for the Ndt ^ Nd s-process. Furthermore, 

("''r^' (-) 

is the propagator of the two intermediate auxiliary fields. 

How to calculate higher-order corrections? The re-summation of all interactions proportional 
to y and A^/j into a Faddeev equation is mandatory since the power-counting classifies 
all diagrams which are built only out of them as contributing equally strong to the final 
amplitude. Ref. proposed to calculate the higher-order corrections by expanding the 
kernel and inhomogeneous part of the integral equation in powers of Q to the desired order 
of accuracy, and iterate then by inserting it into a Faddeev equation {partially re-summed 
N^LO- calculation). This re-summation of some higher-order effects does not increase the 
accuracy of the calculation, which is still set by the accuracy to which the kernel is expanded. 
It does however simplify the numerical treatment, as no divergences in the amplitude or in 
three-body forces Ti. are encountered as the numeric cut-off is removed. It also leads to 
a simple, analytical argument at which orders three-body forces with derivatives enter, as 
discussed below. In contradistinction, a strict perturbation around the LO solution soon 
becomes cumbersome numerically, as full off-shell amplitudes need to be computed and the 
numerical integrals soon start to diverge, making a numerical renormalisation necessary [2^1 • 

The computational effort to solve the integral equations ()3.1f13.4|) numerically is triv- 
ial, as all potentials are separable. A simple Mathematica-code can be down-loaded from 
http://www.physik.tu-muenchen.de/~hgrie. Here, a step-function cut-off A in momen- 
tum space was chosen. The cut-off dependence of the results is discussed below. To reduce 
the numerical instabilities from the poles of the two-nucleon amplitudes and logarithmic 
singularities of the projected nucleon propagator in the kernel, the integral equations are 
- following Hetherington and Schick j2Zj - first solved on a contour in the complex plane. 
The amplitudes on the real axis are then re-constructed by another use of the equations. A 
grid of 70 points does more than suffice for numerical stability. For example, the imaginary 
parts of kcot6 vanish below the deuteron-dis-integration threshold easily to 1 part in 10^. 

Finally, the scattering phase-shift of the Ith partial wave in the quartet and doublet 
channel is related to the renormalised on-shell amplitudes by 

Mfccot5f-ifc' ''^^ M kcot6^l^-ik ' ^ ' 

where x,y = s,t label the matrix entries in cluster-configuration space, and 

fi" = 25> with 2 := J^) (3.8) 
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is the renormalised doublet-amplitude and its wave-function renormalisation. In the doublet 
channel, the only observable process is nucleon-deuteron scattering, Ndt Ndt, i.e. x = 
y = t. Results are not discussed for the un-physical processes Ndt ^ Ndg and Ndg — * Ndg. 
The elastic differential cross-section from A^rf-scattering in the cm-frame is finally EF 



da 



P;(cOSI 



1=0 



k cot 6 



(0 

dM 



ik 



Pi{cosi 



1=0 



k cot d, 



ik 



(3.9) 



A maximum angular momentum 



4 serves to converge the differential cross-section. 



The recipe to compute A^rf-scattering up to N^LO is hence as follows: Expand the auxiliary- 
field propagators and their wave-function renormalisations in the integral equations 1)3.11 / 13 .411 
to the desired level of accuracy in Z-parameterisation as in ()2.19|/IT^m) . and then solve the 
Faddeev equation. To N^LO, the projected nucleon propagator /C-'^ is unchanged. 

Nucleon-deuteron scattering is to N^LO thus completely determined by four simple ob- 
servables of A^A^-scattering: the deuteron binding energy and residue (or effective range), 
and the scattering length and effective range of the ^So-channel. Only the ^Si-channel has 
further unknowns, namely the strength of the three-body interaction Hq at LO and NLO, 
and in addition of H2 at N^LO, as discussed momentarily. They are determined by its 
measured scattering length ad [221 and the triton binding energy Bd, respectively: 

aa = (0.65 ± 0.04) fm , 3^ = 8.48 MeV (3.10) 

But why are two and only two three-body force needed at N^LO? 



3.2 Three-body Force and Lagrangean 

We now turn to the three-body force Ti. already included in the doublet-channel Faddeev 
equation ()3.4|) . The Pauli principle in the quartet channel and the centrifugal barrier in all 
partial waves / > forbids three-body forces in nearly all partial waves at N^LO; see also 
PU] . However, neither rules out a three-body force without derivatives in the ^Si- wave, 
the physically most interesting channel, which contains the triton and ^He as real bound 
states. While one would naively have guessed that such an energy-independent three-body 
force scales as Q° and hence enters only at N^LO, an unusual renormalisation of the three- 
nucleon system in the triton channel mandates its inclusion into the Faddeev equation as a 
LO term As was first pointed out by Minlos and Faddeev [21], this channel suffers 

from a peculiar cut-off sensitivity of the on-shell amplitudes if three-body forces are absent. 
The kernel of ()3.4|) is for / = at short distances identical to the one of an attractive 1/r^- 
potential between the auxiliary fields and the nucleon ^U]. Without a three-body force, the 
wave function would hence collapse, and all observables would become sensitive to Physics 
at very short distances, a phenomenon well-known as Thomas effect Its mathematical 
origin lies in the fact that in the absence of a three-body force, the solution to the integral 
equation is not unique because it allows for a zero mode, as the kernel of the doublet-S wave 
Faddeev equation is not compact [33] . 
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In praxi, the integral equation is solved numerically by imposing a cut-off A, which should 
not be confused with the breakdown-scale of EFT(^). In that case, a unique solution 
exists in the ^Si -channel for each A and 7i = 0, but no unique limit as A ^ cxd. The Nd- 
scattering length, for example, can be tuned as a function of A to have any value between 
minus infinity and plus infinity. As long-distance phenomena must however be insensitive 
to details of the short-distance Physics (and in particular of the regulator chosen), Bedaque 
et al. O CHI m] showed that the system must be stabilised by a three-body force 



^(^'^) = A^E ^2n(A) 



n=0 



It 



A2 



2/7o(A) , 2i72(A) 



A2 



A4 



{ME + 7^ 



(3.11) 



which absorbs all dependence on the cut-off as A ^ oo. It is analytical in E and can be 
obtained from a three-body Lagrangean, employing a three-nucleon auxiliary field analogous 
to the treatment of the two-nucleon channels, see ^1] and the end of this Sub-Section. 

H^n is dimension-less but depends on the cut-off A in a non-trivial way, as a renormali- 
sation-group analysis reveals: Instead of approaching a fixed-point as A — ^ oo, it shows an 
oscillatory behaviour known as "limit cycle" ^21 IMl IHI • 

As one needs a three-body force at LO, Hq ~ to prevent the system from collapse, 
all three-body forces obtained by expanding 7i in powers of E are also enhanced, with the 
interactions proportional to entering at N^^LO fl]. Since a numerical verification of this 
analytical observation is one of the prime advantages Z-parameterisation has, the argument 
is recalled here: Bedaque et al. jJSl noted that a re-parameterisation of the amplitudes 
relates the problem for the LO three-body force to a particularly interesting symmetry, and 
this was extended by Ref. jHj. Building the linear combination 



t 



Wigner 



(0) 

Wigner, 



(0) 

Wigner, 



1 A -A 

m 

■d,ts 



2 VI 1 



(3.12) 



the Faddeev equation becomes in the ^Si -channel 



}C^^\E;k,p) ( \ )+n{E;A)(l 



(3.13) 



oo 

2 f , f2IC(^){E;q,p)+n{E;A) 



<1,P) 



X 



S(E-^,g) A{E-^,q) \ ^-^0) 



%isner(-^; ^' l) 



A 



i(A + Ds) = ^tiV is the ' 'average" A^A^-S-wave scattering amplitude, 
= |strP parameterises the degree to which the ^Si- and ^So-channel differ. 



While E : 

l(A - D 

In the Wigner-S'f/(4)-limit of arbitrary combined spin- and iso-spin rotation [20], the two 
amplitudes are identical, ■jt = Is, Pnt = Pns Vn. In that case, A = 0, and the two equations 
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for ivv"igner,=F decouplc, as first observed in [THj. The essential observation is now that for 
the UV-behaviour of the amphtude, in which all scales are discarded except for the off- 
shell momentum g ~ A in the loop, the Wigner-S'f/ (4)-limit is recovered automatically. 
^Wigner+ obcys the same integral equation as the quartet-S-wave in (|3.ip . where the Pauli 

principle forbids three-body forces without derivatives. Only ^wigner ,- subject to the 
Wigner-5'[/(4) symmetric three-body force Ti. Its integral equation is the same as for three 
spin-less bosons, whose wave-function is well-known to collapse in the absence of three-body 
forces [T21ISI1E31- Hq therefore must be LO to ensure that the physical on-shell amplitude 
is cut-off independent. 

To determine the running of one expands the quantities of ()3.13p for g A;,]?, 7, . . . : 



LO NLO N^LO 

4 1 47 Z -I AME + 872 flZ - 1 V3(Z - 1 



^2 



2g2 



.n^. ^ 1 ME 

IC^'KE;q,p) - ^ + ^4 ' +••• 

^7^(F 2i7o^O(A) , 2<LO(A) ^ 2Hf^^{A) , ME + ^' 2H,{A) , 

^nE;A) - + + + A^ + --- 

(3.14) 

At LO, Ho is independent of two-body observables and is fixed by a three-body datum. 
At NLO, only S has non- vanishing contributions. They depend only on the low-energy 
two-body observables 7 ~ Q and po ~ Q^, but not on the total cm-energy ~ - or 
equivalently, not on the on-shell momentum k Q. The cut-off dependence they induce can 
hence be absorbed by H^^^{K), a momentum-independent correction to the LO three-body 
force whose value now depends also on the two-body observables 7 and po- 

At N^LO, more corrections arise which depend on 7 and po. In addition, both the 
projected nucleon propagator K, and the iViV-scattering amplitude S show corrections pro- 
portional to ME. To match the behaviour of the off-shell (UV-) amplitude on the physical 
on-shell momentum E., the inclusion of an energy- dependent three-body force H2 is thus 
mandatory. Its value is not determined by two-body observables. The expansion in powers 
oi E proceeds straightforwardly: Another three-body force H2n{A) proportional to 

E"^ is needed every even order in Q from expanding both S and /C'-^-'. if2n(A) first enters 
at N^"LO, independent of two-body observables. To determine its strength, one additional 
three-body datum is needed. The power-counting for the three-body forces is hence 

Ho ~ , H2 ~ Q-^ , H2n ~ . (3.15) 

Corrections to the limit q ^ 00 from the breaking of Wigner-5'[/(4)-symmetry can only 
induce a momentum-dependent three-body force, because the only three-body force without 
derivatives is necessarily Wigner-S'f/(4)-symmetric fH]- They can be shown to be of higher 
order, see HH for a more detailed discussion. 
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The chain of isotropic (S-wave) Wigner-S'?7(4)-symmetric three-body forces is also found 
from a three-nucleon auxihary-field Lagrangean, in complete analogy to the treatment of 
the two-nucleon auxiliary fields in Sects. I2.1f12.3l [H] : 



3N 



2/3(A) 



[t^{{a^N)dl-iT^N)df)+E.c.] 



(3.16) 



2 

n=l 



6M M 



where the auxiliary field t has the quantum numbers of the triton/^He, i.e. spin and iso- 
spin i. The relative coupling strength between the processes t Ndt and t — > Nds is 
fixed because there is only one three-body force without derivatives, which also happens to 
be Wigner-S'[/(4)-symmetric [T!^ I14j. This field is treated analogously to the two-nucleon 
auxiliary fields: At LO, the propagator of a "triton" with kinetic energy E and momentum 
p is proportional to ^, with insertions proportional to h2n{^) suppressed by n powers of the 
triton kinetic energy: 



n 



(3.17) 



Comparing the final result of the resultant three-body force in the ^Si-wave ()A.13|) with 
()3.4|) . one finds that the couplings 1/3 (A) and /i2n(A) are in the cm- frame related to Tt by 



n{E-K) 



2/2 Mn 



E 



m,n=l 



'-^^-1(ME + 7,Y' 



(3.18) 



or matching order by order the dimension- less three-body strengthes H2n{^) of ()3.1H) . 

^„,A) = -IM ^ . „.M = -IM ^ UM) ..... (3.19) 



2/2 Mfi ' ^ ' ?/2 ^MQ^^2 -V 

SO that the couplings are particularly simple when one chooses 



(3.20) 



as they follow then from ()3.15|) the natural scaling laws - with the exception of 2/3 (A): 

2/|(A) ~ , 1] ~ Q° , /i2n(A) ~ g° (3.21) 

Notice that in contradistinction to in Sect. 12.11 the power counting does not require 
even a partial re-summation of the triton propagator. 
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3.3 The 2Si- (Triton-) Channel 



That Z-parameterisation improves convergence is most importantly seen in the only channel 
in which a bound state is found, and in which the three-body force enters already at LO to 
stabilise the system from collapse. In particular, it can be instrumentalised to confirm that 
the first momentum-dependent three-body force enters at N^LO, as the analytical argument 
of the preceding section demands. 

Previously, Hammer and Mehen calculated the ^Si-wave phase shifts below the deuteron 
break-up to NLO by sandwiching the ERE- correct ions to the two-particle propagators be- 
tween the LO half off-shell amplitudes [20] • ERE-parameterisation was also used in the 
recent article which proposed partially re-summed N"LO-calculations to simplify numerics 
and showed that the momentum-dependent three-body force H2 is conceptually necessary at 
N^LO J3]. An estimate of NLO-effects in Z-parameterisation by Afnan and Phillips jSZI in- 
cluded the part of the two-particle propagator (j2.19|/IT^m) which leads to the correct residues 
Zs/ti but not the part which changes the off-shell propagation of the deuteron field, iterating 
again the perturbatively expanded kernel. 

In the following, one has to differentiate between two kinds of convergence: First, the 

results should of course agree with available measurements, at the level of accuracy predicted 

by EFT. Second, it is vital for a reliable error-estimate and hence for predicting the accuracy 

of the calculation that EFT can demonstrate that contributions which are classified as 

higher-order are indeed suppressed, i.e. that corrections from order to order become smaller 

in the range of validity. In particular whenever no or scarce data are available, this a-priori 

error-estimate is essential to minimise the theoretical bias in the description of few-body 

properties. Z-parameterisation improves both variants over ERE-parameterisation. 

The ERE- and Z-parameterisation results for the phase-shifts of the ^Si-channel at LO, 

2 ,__ 

NLO and N^LO are shown in Fig. El together with the only available phase-shift analysis [SHj . 
dating from 1967, and results using the Argonne V18 potential supplemented with the 
Urbana IX three-body force l4Uj . The bands in the N^LO-results come from varying 
the cut-off in the Faddeev equation between 200 and 900 MeV, when A-independence is 
achieved within drawing-accuracy, except at those points at which the strengthes of the 
three-body forces diverge, see later Fig. El 

At this order, considerably more sophisticated and involved potential-model calculations 
must agree with the EFT-predictions, if all are fitted to the same low-energy observables. 
For lack of a direct phase-shift analysis of the three-nucleon system, convergence to experi- 
ment is discussed by comparison to such a potential-model calculation which reproduces the 
same two-nucleon observables, triton binding energy and scattering length in the ^Si-wave. 
At higher orders, EFT-calculations may be more accurate than those of potential models, 
because the EFT-result will depend on new three-body low-energy coefficients which do not 
enter in potential models, while a systematic treatment of three-body forces is built into 
EFT(7f). 

At LO, both parameterisations are identical. As expected from the discussion in Sect. 12.21 
the NLO corrections are larger for Z-parameterisation than for the ERE-version due to 
the inclusion of the full strength of the deuteron residue. Higher-order corrections in Z- 
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Figure 3: Comparison of the real (top) and imaginary (bottom) parts of the neutron- 
deuteron ^Si -phase-shift at LO (dotted), NLO (dashed) and N^LO (dark band) in ERE- 
parameterisation (left) and Z-parameterisation (right) as function of the cm-momentum 
k. The band in the N^LO -curve shows the variation induced by shifting the cut-off from 
A = 200 MeV to oo. Dots: phase-shift analysis from 1967 jHHj; crosses: Argonne V18 results 
including the Urbana IX three-body force from jSH] below, and from jlU] above break-up. 
At LO and NLO, the results for A = 900 MeV are shown. The dot-dashed curve in the 
imaginary parts of the ERE- result is a partial N^LO-calculation as described in the text. 



parameterisation should be considerably smaller than in ER15-parameterisation, because 
their typical scale from A^A^-scattering is in the former set by {Zt — l)/2 ^ 0.3, while 
in the latter, it is 'jtPot ~ 0.4. This amounts to a difference in the predicted accuracy 
of the N^LO-calculation of (0.3)^ ~ 3% in Z-parameterisation versus (0.4)^ ^ 6.5% in 
ERE-parameterisation for typical momenta in the three-body problem. Indeed, the cor- 
rection from NLO to N^LO is in general larger for the ERE-parameterisation than for 
the Z-parameterisation, and the result closer to the potential-model values. The plot of 
the logarithmic deviation of the on-shell point of the inverse /T-matrix, k cot 6, from the 
potential-model numbers at each order in both parameterisations verifies these findings on 
the quantitative level, see left panel in Fig. |3] 
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In particular, the ERE-result for the imaginary part of the phase-shift has at large 
momenta not converged to the potential-model results at N^LO, while Z-parameterisation 
agrees with these computations. As the effective-range effects are expected to be the dom- 
inant corrections in ERE-parameterisation, the imaginary part of the ERE-result in Fig. 121 
contains also a computation in which these are included in the deuteron propagator ()2.4|1 
to order {'~ftPotY- This is only a partial N^LO-calculation, neglecting shape-parameter, rela- 
tivistic, SD-mixing and other corrections, but improves the agreement to the level achieved 
already at N^LO in Z-parameterisation. Close to k = ^ A^, both parameterisations 
become un- reliable as the NLO and N^LO-corrections are comparable in size. 




30 50 70 100 10 15 20 30 50 70 100 

k [MeV] k [MeV] 



Figure 4: Convergence of k cot 6 as function of the cm momentum. Left: deviation from the 
result of the potential-model calculation [211 EHl- Right: dependence on cut-off variations 
between A = 200 and 900 MeV. Dotted: LO; dashed: NLO; sohd: N^LO; thin lines: 
ERE-parameterisation; thick lines: Z-parameterisation. Thick dots on the right: "N^LO"- 
calculation in Z-parameterisation with H2 set to zero. 



At least as important is that the internal convergence is drastically improved in Z-parameter- 
isation. This is evident from the right side of Fig.Hlwhich compares the uncertainties induced 
into the N^LO phase-shifts by cut-off variations. As outlined in the Introduction, the main 
strength of the EFT-approach is that it allows for an a-priori estimate of the theoretical 
uncertainties of a calculation because contributions to the amplitudes are ordered by a small 
expansion parameter Q. N^LO corrections to kcot6 should typically be of the order 

A(A;cot<5) ~Q"= (^^^ (3.22) 

compared to the LO result. Typical low-momentum scales ptyp in the three-body system are 
the binding momenta of the two-nucleon real and virtual bound states, 7s ~ —8.0 MeV, jt ~ 
45 MeV and the scattering momentum k. In addition, the three-body forces are determined 
in part by the typical three-nucleon bound-state momentum 7^ ~ a/M5^ ^ 90 MeV, 
the triton binding energy. The breakdown scale A^ ^ of the theory is the scale at which 
higher-order corrections become comparable in size. Like the actual size of the expansion 
parameter Q, its value must be verified in actual calculations. 
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As observables at low energies must be independent of details of short-distance physics, 
they must be independent of the arbitrary regulator A up to the order of the expansion. In 
other words, the physical scattering amplitude (or more accurately its physically relevant 
part, kcotS) must be dominated by integrations over off-shell momenta q in the integral 
equations ()3.1|/IH.4|1 in the region in which the EFT is applicable, q < A^. As argued e.g. by 
Lepage [l^, one can therefore estimate sensitivity to short-distance Physics, and hence 
provide a reasonable error analysis, by varying the cut-off A between the break-down scale 
Ajf and oo. In praxi, an upper limit A^ax = 900 MeV suffices since the phase-shifts are 
essentially cut-off independent beyond A ^ 600 MeV. 

On the right in Fig. 01 the logarithmic cut-off variation is displayed as function of the 
cm momentum, both for ERE- and Z-parameterisation. Several points are worth noticing. 

First, the cut-off variation is substantially smaller in Z-parameterisation over the whole 
range. For example, at the typical three-body scale k ~ y/MB^ ~ 90 MeV, the cut-off 
variation is decreased from 15% and 4% at NLO and N^LO in ERE-parameterisation to 
12% and 1.3% in Z-parameterisation ^. 

With k the dominant low-energy scale for momenta above the breakup point (47^^/3)^/^ ~ 
52 MeV, one can secondly verify the estimate (|3.22p for the corrections by fitting n to the 
nearly straight lines in the right panel of Fig. |3] in the momentum range between 70 and 
100 to 130 MeV. At lower momenta, the slope does not change significantly from order to 
order because the dominant corrections are ~ (7s^f/A)" and hence nearly independent of k. 
The power-law obtained for each order and parameterisation is listed in Table ^ The slope 
increases in Z-parameterisation by one unit from LO to NLO, and by two units from NLO to 
N^LO. The latter may stem from the partial inclusion of higher-order graphs in the partially 
re-summed integral equation. The ERE-parameterisation follows a weaker power-law than 
Z-parameterisation, indicating again greater cut-off sensitivity and less accuracy. 

Thirdly, extrapolating the relative errors from the fit region, the cut-off variations at NLO 
and N^LO become comparable at /c ~ 150 — 200 MeV, matching the estimate A^ ~ 771.,^. 



order 


n (Z-param.) 


n (ERE-param.) 


LO 


1.8 - 2.1 


NLO 


2.9 


2.6 


N^LO 


4.8 


3.7 


"N^LO" , if2 = 


3.1 


2.8 



Table 1: Fit of the logarithmic cut-off variation 1 — nnn^An of A; cot 5, shown in 

° fccotd(A=900 MeV) ' 

Fig. m for momenta higher than 70 and lower than 100 to 130 MeV to a power-law (fc/A^)"' 
for ERE- and Z-parameterisation at LO, NLO and N^LO, and at N^LO without momentum- 
dependent three-body force H2. Varying the upper limit changes only the numbers for LO. 



''^The deviation at N^LO to the potential-model results is at the same scale improved from 9.1% to 5.7%. 
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Finally, one confirms numerically the analytical finding in ^1], recalled in Sect. 13 .ll that 
a momentum-dependent three-body force H2 enters at N^LO to increase accuracy. Setting 
H2 = does not only substantially increase the cut-off dependence of the N^LO-calculation, 
putting it close to the NLO dependence. As the slope at large momenta is also practically 
unchanged compared to the NLO value, H2 is necessary to improve cut-off independence. 
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Figure 5: Dependence of the dimension-less three-body couplings at N^LO on the sharp 
momentum cut-off A G [140; 2800] MeV in Z-parameterisation. Triangles: Hq{A); crosses: 

i^2(A). 



Notice that observables are at N^LO converged to be practically cut-off independent as 
low as 600 MeV, while the dimension-less three-body forces Hq and H2 vary dramatically 
even in a small window of cut-off variations, see Fig IHl Only when the absolute value of 
Hq or H2 is unnaturally large (> 10) is cut-off sensitivity felt. In that case, more of the 
string of three-nucleon forces would be needed. Such singular points (e.g. atA ~ lOOOMeV) 
are not considered here. The actual numerical values of Hq and H2 are quite sensitive to 
the regulator chosen, and to the precise values for 'jg/t and Zg/t used. It is thus not the 
three-body forces which are large but the effect of variations in the cut-off A on observables. 
No simple connection between Hq and H2 exists at N^LO, except that they both show a 
limit-cycle with the same period, which depends on 7^/^ and Zg/f 

To complete the discussion. Fig. El displays the NLO-prediction for the Phillips line 
in Z-parameterisation, compared to various old and modern potential model calculations 
which share the same two-body on-shell Physics. In the three-body system, their widely 
different off-shell behaviour predicts however different three-body bound state energies and 
Nd scattering lengthes. The correlation between scattering length and bound state energy 
was first discussed by Phillips jHj and explained by Efimov |12]- The spread induced by 
varying the cut-off between 140 and 900 MeV covers the off-shell dependence of all potential 
models. Note the pole in the scattering length at B'^ ^ 36 MeV at NLO for a cut-off 
A = 900 MeV. This intrusion of another bound triton state is however outside the range 
of applicability of EFT(/). As hyp ^ 180 MeV > A^ and the widthes of the LO and NLO 
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Figure 6: Prediction for the Phillips line at LO (dotted; lightly shaded region) and NLO 
(dashed; dark shaded region) in Z-parameterisation when the cut-off is varied in the region 
A e [140; 900] MeV. Left: The dots correspond to the predictions for the triton binding 
energy and doublet scattering length in different models with the same two-body scattering 
lengths and effective ranges as inputs jl21- experimental result. The outcome in various 
modern high-precision A^A^-potentials is indicated by crosses jlSl- Right: Larger range, the 
rectangle showing the region covered on the left. 



band are comparable, the EFT-calculation has not converged there. 
3.4 Including the Effective Range to All Orders? 

Why not re-sum all effective-range parameters up to a given order and replace the expanded 
version of the two-body scattering amplitudes (|2.12j) with the version ()2.6f12.2l|) . 

Ds/tipo,P) 5S ^ , , 

Is/t - '-^{1% + f - Mpo) - E PnslAlsIt + f - Mpo)"+l - Vf - Mpo 

n=l * 

(3.23) 

possibly with all coefficients after Pnsjt set to zero at N^'^+^LO? This would immediately 
restore all effective-range corrections, give the correct pole residues, and is mandatory if 
Po7 ~ 1- With this re-summation, the scattering length in the ^Ss-channel and phase-shifts 
of the higher partial waves [IHl CHI EOl 122 ES] were found with an accuracy of < 4%. It was 
later also shown to simplify analytical calculations in the two-nucleon system |21j . 

However, the improvement in accuracy is only superficial: The higher ERE-parameters 
Pnsjt are known with less and less accuracy, and other higher-order terms like P-wave inter- 
actions and relativistic corrections are still not included. 

More severely, the propagator (|3.23|) has additional spurious poles, see Table |21 one for 
each S-wave at NLO (pos/t 7^ 0), and two more at N^LO (pi^/t 7^ 0). Their positions change 
as more parameters Pnsjt are included, and their residues add with the residue of the low- 
lying real and virtual bound state to zero. Albeit they all lie strictly speaking outside the 
range of validity of EFT(]f), 7spur,s/t > A^f, they are in particular in the ^So-channel close 
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enough to influence observables at momenta k above the break-up threshold, in particular 
as the closest of the spurious bound states have rather large residues. Recall that the 
half-off-shell amplitude goes to zero in the ^Si -channel only slowly, A{k,p) ocl/p [121 Ell- 





■^Si-channel 


^So-channel 


order 


7spur,i [MeV] 


•^spur,i 


7spur,s [MeV] 


■^spur.s 


NLO 


178 


-1.69 


152 


-0.90 


N^LO 


-385 
170 ± 108 i 


-0.14 
-0.77^0.60 i 


133 
-62 ±351 i 


-0.71 
-0.10^0.12 i 



Table 2: Binding momentum 7spur = \/ MBspur and dimension- less residue analogous to 
()2.10|/IT^n|l of the additional spurious bound states in the ^Si- and ^So-scattering amplitudes 
when the effective-range parameters are re-summed as in ()3.23|1 . 

In the Nd-sjstem, these additional poles generate new cuts in the solution to the integral 
equations: In the doublet channels, they conspire to pose an additional technical problem 
when one performs a rotation of the integration contour into the complex plane in the 
Hetherington-Schick procedure. On the other hand, expanding the two-particle propagator 
to a given order as in ()2.19|) . only the physical pole with •js/t ^ is found in the two-particle 
scattering amplitude, and all problems with spurious bound-states are easily avoided. 

Most severely, the asymptotic behaviour of the two-nucleon propagator is changed from 
1/p to at large p ^ 1/pos/t- This leads to a non-trivial modification of the argument 
for the inclusion of a three-body force in the ^Si-channel, as the resulting integral equa- 
tion is then not analytically soluble in the UV-limit. Therefore, it is not straight-forward 
to construct analytically the cut-off dependence of the three-body force as in [T^ IT^ and 
demonstrate that it solves the Thomas problem of an attractive 1/r^-potential. Gabbiani 
argued that with the kernel becoming less attractive, behaving only like 1/p^ (i.e. in po- 
sition space like the Coulomb potential 1/r) , the cut-off sensitivity of the Wigner-S'f/ (4) 
symmetric component ()3.13|) of the integral equation in this channel disappears jlHl- Thus, 
he claims that no limit cycle is found, and no three-body force is needed to achieve cut-off 
independence. According to him, the resulting phase-shifts without three-body forces in the 
^Si -channel fail however miserably. I cannot confirm these findings but see that the phase- 
shifts in the ^Si -channel are markedly cut-off sensitive even at high A, and that three-body 
forces are necessary when the effective range is included to all orders. These findings will 
be discussed in a sub-sequent article jlTj. 

Re-summing the effective-range corrections poses therefore numerically as well as techni- 
cally non-trivial problems most importantly in the i-wave (but also in the higher doublet- 
channels) while no increase in accuracy is achieved. To insert on the other hand the expanded 
two-nucleon amplitudes ()2.19|) in the kernel of the integral equation, as proposed in jj^ and 
done here is a straightforward way to calculate phase-shifts to high accuracy, free of the 
problems mentioned above. In the end, all regularisation and re-summation schemes - in- 
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eluding those whieh generate spurious bound states outside the range of vahdity of the EFT 
- must agree to the order of accuracy of the calculation. They must also have the same num- 
ber of independently fixed three-body forces. To set up a consistent power-counting scheme 
and computations with good numerical convergence in a particular version can however be 
cumbersome and time-consuming. 

3.5 More Partial Waves 

After discussing the merits of Z-parameterisation in the only channel in which a three-body 
force enters already at LO and which exhibits a bound state of three nucleons, we now turn 
to the other partial waves. The solutions of the integral equations ()H.1|/I3.4|1 for the real and 
imaginary parts of the phase-shifts are displayed in Fig. [7| for the doublet channel, and in 
Fig. El for the quartet channel. The quartet-S wave is shown separately in Fig. |H1 In the 
following, they are compared to potential-model calculations, and their convergence from 
order to order is discussed and juxtaposed with the phase-shifts obtained when the effective 
range is included to all orders. Table El lists low-energy parameters of lower partial waves. 

Previous work on the higher phase-shifts in the A'^d-system in EFT(^) focused mainly 
on the '^S 3 -channel, where prf-scattering jl^ and the scattering length ^Sl ^1 120] and 
phase shifts [121 (201 122] of nd-scattering were computed up to N^LO. The phase-shifts of 
the other partial waves were presented in Ref. [221 to N^LO. The NLO-results of [201 and 
[22] were obtained in "strict perturbation", i.e. sandwiching the NLO-correction to the 
deuteron propagator proportional to CQs/t once between the LO-solution to the half-off- shell 
amplitude. Ref. [221 used already the Z-parameterisation to fix cot at NLO. In each case, 
the N^LO-calculations were carried out with the re-summed deuteron propagator (|3.23p . 

Here, we use as in the ^Si -channel Z-parameterisation at each order and include per- 
turbatively all corrections up to NLO and N^LO respectively into the kernel of the integral 
equation. As already discussed in Sect. 13. ll this contains some higher-order effects and leads 
to technical simplifications but - like a re-summed deuteron propagator - not to increased 
accuracy. The goal is again to show that convergence and agreement with results from 
potential-model calculations is improved. 

Like in the discussion of the ^Si-wave in Sect. 13.31 the N^LO-results are compared to 

2 

phase-shifts of a modern high-precision A^A^-potential. Numbers for nrf-scattering using the 
AV18-potential with the Urbana IX-three-body- force were reported by Kievsky et al. [S^ EO] 
and Hiiber et al. [IHllHl- Within the drawing accuracy, they agree with the phase-shifts 
reported for the Nijmegen potentials 93, I and II, and for CDBonn [lU]- As the partial waves 
do not mix in EFT at N^LO, the potential-model results for the ^^"""^/j -partial waves were 
grouped into bins of fixed spin s and angular momentum /, with total angular momenta 
J = Z ± I in the doublet channel, and j = Z ± {i; |} in the quartet. 

The higher the partial wave, the better the agreement with the potential model calcu- 
lations. This comes as no surprise, since the strong centrifugal barrier for large angular 
momenta eliminates sensitivity on short- distance physics and the solution to the integral 
equation approaches the result of the Born approximation. 
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Figure 7: Real (left) and imaginary (right) part of the phase-shifts of the higher partial 
waves in the doublet channel in Z-parameterisation as function of the cm momentum k, 
A = 1000 MeV. Dotted: LO; dashed: NLO; solid: N^LO; dot-dashed: effective-range 
corrections summed to all orders. Result of AV18 + U IX from ^\ below break-up, from j3H] 
above, and at /c = 104 MeV from ^Ij wave denoted by x: j = / — |; +: j = I + ^. 
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Figure 8: Phase-shift in Z-parameterisation (left) and comparison of ERE- and Z- 
parameterisation for kcot6 (right) of the ^Ss-wave as function of the cm momentum k; 
A = 1000 MeV. Notation as in Figs. [7|and|ll and in addition on the right: boxes/stars: 
partial wave analysis [SHI (for points with error-bars as reported in jS]); dot at A; = 0: 
measured scattering length [221 • 



quartet-S 


a [fm] 


ro [fm] 


order 


Z-param. 


ERE-param. 


Z-param. 


LO 


5.091 


-0.1 


NLO 


6.410 


5.938 


1.95 


N^LO 


6.354 ±0.020 


6.204 


1.8 ±0.1 


experiment [20] 


6.35 ±0.02 




modern A^A^-potentials [43, 50^ 


6.34... 6.35 




N^LO re-summed, cf. [H Ci 120] 


6.365 


1.8 





doublet-P 


doublet-D 


quartet-P 


quartet-D 


order 


kq [fm^] 


Ki [fm ^] 


kq [fm^] 


Ko [fm^] 




kq [fm^] 


LO 


33.1 


-3.2 


-268 


-76.7 


1.0 


531 


NLO 


53.9 


-2.1 


-456 


-139 


0.45 


894 


N^LO 


53.3 ±0.2 


-2.1 


-456 


-140 


0.43 ±0.01 


894 


N^LO re- summed 








-140 


0.43 


894 



Table 3: The effective-range parameters ()3.24|) of the low partial waves in the A^rf-system 
from Z-parameterisation in the first three orders; theoretical accuracy discussed in the text. 
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Figure 9: Real (left) and imaginary (right) part of the phase-shifts of the higher partial 
waves in the quartet channel in Z-parameterisation as function of the cm momentum k, 
A = 1000 MeV. Notation as in Fig. d and □: j = / - |; 0: j = / + | for the '^Ij-wave. 
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The EFT-calculation reproduces of course the well-known behaviour of the phase shifts at 
small momenta, e.g. |28j . 



1 (0 

^2^+1 cot ^(0 ^ + + 0{k^) VA; ^ , (3.24) 

4^ 2 

where the S-wave parameters Kq^^ and nf'^ are the scattering length a and effective range 
ro, the P-wave parameter Kq^"* the scattering volume, etc. Table El lists parameters for low 
partial waves of the A^rf-system. The theoretical accuracy may for example be estimated 
conservatively by Q ~ ^ ~ | of the difference between the NLO- and N^LO-result, or by 
the difference to the N^LO-result with re-summed effective-range corrections as in (|3.23p . 
since all N^LO-calculations must agree within the N^LO-accuracy. The numerical uncertain- 
ties are easily covered by these error-bars. In the doublet channel, no re-summed results are 
given due to the non-trivial cut-structure induced by the spurious poles of the two-nucleon 
propagators, see Sect. 13.41 Furthermore, no numbers are presented for the ^Si -channel as 

2 

the measured scattering length ()3.1()j] is used to determine the strength of the three-body 
force, and the effective range is very large, ~ 500 fm, rendering an effective-range expansion 
useless for this partial wave. In the other partial waves, the parameters involve the typical 
length scale ~ 2 ... 5 fm, so that an effective-range expansion is often useful for momenta 
up to the deuteron break-up, A;max ^ 20 ... 52 MeV. More effective-range parameters are 
obtained without difficulty, but they are less likely to be measurable in the upcoming partial 
wave analysis of the nd-system |51j . 

The quartet-S wave scattering length has drawn substantial interest recently as its knowl- 
edge sets at present the experimental uncertainty in a recent, indirect determination of the 
doublet scattering length ()3.10p |52j . At N^LO, the result of Z-parameterisation in EFT{f) 
(using a perturbatively expanded kernel which is the inserted into the integral equation), 
[6.35 ±0.02] fm, agrees very well with experiment [221 and modern high-precision potential- 
model calculations (which give = [6.34. . .6.35] fm) 150j. albeit partial-wave mixing, 
iso-spin breaking and electro-magnetic effects are not present in EFT(^) at N^LO. is 
obviously to a very high degree sensitive only to the correct asymptotic tail of the deuteron 
wave function, as comparison between the NLO result with its correct deuteron residue, the 
N^LO-result and the result of the re-summed deuteron propagator shows. As the amplitude 
decays at large off-shell momenta as i/|)3.i7... this is not surprising. 

Comparing k cot 6 in this partial wave at low energies in Fig. |S[ the ERE- version con- 
verges much slower than the result of Z-parameterisation. For example, the N^LO-scattering 
length differs from experiment and the Z-parameterisation result of the same order still by 
2%, and the correction from NLO to N^LO is with 5% much larger than the correction in 
Z-parametrisation. Indeed, partial N^LO- and N'^LO-calculations in ERE-parameterisation 
which include only effective-range effects bring one slowly closer to the experimental number: 

a(^S)ERE = (5.09 + 0.84+ 0.27+ 0.10 + 0.04 + . . . ) fm ^ 6.34 fm (3.25) 

LO NLO N^LO "N^LO " "N^LO " 
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The re-summed effective-range result of EFT(/) at N^LO was first reported in [THl ll9j. The 
ERE-result at NLO in [20], 6.7 fm, was found not by re-summing the expanded kernel as 
here, but by inserting the NLO-correction only once between the LO amplitudes. Thus, 
convergence is also sped up by iterating the perturbed kernel. 

EFT shares also the rather complex structure found even at low k, e.g. in the doublet- 
and quartet-D and -F waves. Some disagreement with the potential-model results is seen 
in some imaginary parts at rather high momenta k > 120 MeV, close to the breakdown 
scale of EFT(^), where the expansions parameter Q is approaching 1. Given that the 
imaginary parts are very small compared to the real parts and also most sensitive to details 
of the break-up reaction nd — > nnp, the overall agreement is not dis-encouraging. The 
most notable deviation occurs in the doublet-P wave, where internal convergence of the 
EFT-result is achieved up to more that k ^ 100 MeV, while significant deviations from 
potential-model calculations start around breakup. In this as in all other channels, cut-off 
dependence is negligible above A ~ 900 MeV (indicating also that a three-body force is not 
needed for convergence, see [SD]), and numerics is stable. The discrepancy - already found 
in [221 ^ has therefore a physical origin. We may expect that it is cured by cancellation with 
the well-known, rather strong P-wave interaction in the iVA^-system. It enters at N^LO, 
where the splitting between the partial waves will also be reproduced more realistically. 
Clearly, more work is needed here. 

While the NLO-corrections to the LO EFT-result are as expected rather large in Z- 
parameterisation, the N^LO-modifications are in the real parts of the phase-shifts tiny, and 
do not exceed 40% of the NLO corrections in the imaginary parts even as /c — *• ~ 
m^, the ^P-wave being again the only notable exception. The main difference between 
the strictly perturbative NLO-calculation of [22] which already utilised Z-parameterisation 
and the partially re-summed version used here is that now the NLO and N^LO-results lie 
practically on top of each other, while they noticeably differ in strictly perturbative Z- 
parameterisation above break-up. The difference does however never exceed 5%, and hence 
lies within the accuracy of the NLO calculation. The N^LO-results of ^21 1201 122] using 
the fully re-summed deuteron and ^So-propagators ()3.23|) are nearly in-discernible from the 
partially re-summed version ()2.19|) in the real parts even at hight k. In the imaginary 
parts, they do usually not deviate by more than 10% of the correction from NLO to N^LO 
at /c ~ 140 MeV. This lies again well within the power-counting prediction that partial 
inclusion of N^LO-effects will not increase the accuracy of the N^LO-result. It gives however 
a band within which one expects a full N^LO-calculation to lie. 

Finally, an unsolved technical issue has to be mentioned. The problem that the spurious 
poles of the re-summed two-body propagators induce additional cuts in the solution of the 
integral equation in the doublet channel, was already discussed in Sect. 13.41 Such breaches of 
unitarity below breakup, induced by naive integration over a cut at high off-shell momenta, 
are the reason why the phase-shifts of the doublet-P and doublet-D-wave have a small non- 
zero value at the breakup point. As the half-off- shell amplitudes of higher partial waves 
converge substantially faster than of the doublet-S-wave, these violations become however 
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less and less noticeable. For the higher partial waves, they disappear in the numerical noise. 

To summarise, the overall convergence of the phase-shifts is good in all channels, both 
internally and to the available potential-model calculations. The only noticeable exception is 
the doublet-P-wave, in which the radius of convergence is limited to momenta k below break- 
up. Neglecting partial- wave splitting, simple observables like the differential cross-section do 
therefore agree with potential-model results and direct experimental measurements, within 
the level of accuracy of the EFT-prediction. As an illustration, the elastic differential cross- 
section is shown in Fig. for three neutron energies: just below the nd-breakup [k = 
50 MeV), and at two momenta close to the breakdown scale of EFT(/), k ^ 90; 120 MeV. 
The accuracy of the calculation is estimated by varying the N^LO-correction to k cot 6 by 
Q ~ 'jt/'^n ~ 1/3 around its central value. The agreement with experiments |33] 
and potential-model calculations [23 IH] is quite satisfactory even close to the breakdown 
scale k ~ m^f. The main source of deviation is of course the insufficient description of the 
doublet-P-wave. 

More complex variables need of course more refinement. For example, one obtains a 
zero-result for the nucleon-deuteron vector analysing power Ay in EFT(^) at N^LO because 
partial- wave splitting is absent. The well-known under-prediction of this observable in all 
modern, high-precision A^A^-potentials |2] - even when supplemented with three-body forces 
- is a challenge a higher-order EFT-calculation has to meet. 

4 Conclusions and Outlook 

Z-parameterisation ^H] fixes the parameters of "pion-less" Effective Field Theory at very 
low energies up to N^LO, i.e. up to an estimated theoretical accuracy of < 3%, to the 
correct deuteron pole position and residue, instead of the effective-range parameters. As 
the asymptotic fall-off and normalisation of the deuteron wave-function at large distances 
is correctly reproduced, this improves drastically convergence in all processes in which the 
deuteron is either found as in- or out-state, or as sub-cluster of a more complex few-body 
system. The reason for the improvement is found in a rather large expansion parameter of 
ItPot ~ 0.4 for effective-range corrections, so that the deuteron residue is regained only up 
to 7% at N^LO when fitting to the effective range, while it is exact already at NLO in Z- 
parameterisation. Thus, the speed of convergence is faster, while the radius of convergence, 
set by the first scale on which new Physics enters, is of course not changed. 

In this article, Z-parameterisation was extended and applied to deuteron-nucleon scat- 
tering, see Sects. El and 13.11 To perform the calculations, a computationally simple and 
convenient scheme was again employed which first expands the kernel to the desired order 
of accuracy in the power counting and then iterates it to all orders ^3]. 

Numerical problems which origin from additional, spurious poles in the TVTV-amplitudes 
when the effective-range corrections are re-summed to all orders in the deuteron propagator 
are avoided in Z-parameterisation, Sect. 13.41 Therefore, standard techniques are readily ap- 
plicable, and the calculation of scattering-observables also above break-up becomes compu- 
tationally trivial, see a Mathematica-code at http : //www . physik . tu-muenchen . de/~hgrie. 
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e [rad] 

Figure 10: Elastic differential cross-sections of nd-scattering at N^LO in the cm-frame in 
Z-parameterisation (solid) with estimated accuracy (band). Dot-dashed: AV18 + Urbana 
IX from at neutron kinetic energy Eiab = 3 MeV, AV18 from ^ at Eiab = 10; 18 MeV; 
diamonds: rad-data [521 (taken at 3; 10.25; 18 MeV); stars: pci-data 
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In all partial waves, the agreement of the phase-shifts with sophisticated potential-model 
calculations is increased in particular in the real and imaginary parts above the deuteron 
break-up point, see Sect. 13.51 What is more, the corrections from one order to the next 
in the EFT(^)-expansion are smaller in Z-parameterisation, so that a partial re-summation 
of effective-range effects increases indeed the theoretical accuracy of the calculation. This 
allowed also for high-accuracy predictions of effective-range parameters in the low partial 
waves of the nd-system, which can be compared to an upcoming partial-wave analysis of 
the nd-system j^I]. Most notably, the quartet-S wave scattering length is predicted as 
Qq = [6.35 ± 0.02] fm at N^LO, using only four observables from the iViV-system as input, 
namely the asymptotic form and normalisation of the deuteron and virtual ^So-bound state. 
This result compares both in magnitude and uncertainty favourably with the most advanced 
potential-model calculations ([6.34 ... 6.35] fm ED]), and with the experimental value 
([6.35 ±0.02] fm [2ni)- Simple observables like differential cross-sections are well reproduced 
as high as -Eiab ~ 15 MeV. 

The most striking success lies however in the ^Si-channel, which contains the triton 
as real bound state. Sect. 13.31 The numerical analysis of the cut-off dependence of its 
phase-shift in Z-parameterisation clearly supports the analytic finding in Ref. ^3] that one 
and only one momentum-dependent three-body force enters at N^LO, namely the Wigner- 
S'?7(4)-symmetric one with two derivatives. Fixing its strength to the triton binding energy, 
the accuracy of the calculation - deduced from the cut-off dependence of the answer - is 
improved at N^LO from ~ 4% to ~ 1% at a momentum scale of ~ 90 MeV typical for 
the three-body system. The cut-off dependence follows in Z-parameterisation the pattern 
predicted by the power-counting. The difference to potential-model calculations is also 
substantially diminished. 

The major source of uncertainty in the triton channel stems thus not from the theoretical 
accuracy of the EFT(^)-calculation, but from the error in the determination of its nd- 
scattering length, which is known only to ~ 7% accuracy. To reduce this uncomfortably 
large error bar to an accuracy of ~ 0.7% is the goal of an ongoing direct measurement of the 
incoherent scattering length at the Paul-Scherer-Institute jSE]- A recent measurement of the 
bound coherent scattering length |H2] relied for the extraction of ad = [0.645 ±0.003(exp) ± 
0.007(theor)] fm on theoretical input for the ^S 3 -wave scattering length [3Uj . 

Future work includes a complete N^LO-calculation, including iso-spin breaking effects and 
partial- wave splitting and mixing (which should improve in particular the doublet-P waves), 
the extension to include Coulomb interactions for the ^He- and pti-system, and the coupling 
of electro-weak probes to the three-nucleon system. Not only will this allow for direct 
comparison with a cornucopia of data and help to shed light on long-standing puzzles like 
the Ay-problem. It will also lead to predictions with an accuracy relevant for nuclear astro- 
physics and neutrino physics at very low energies, contributing to the ongoing efforts to 
improve our knowledge about big-bang nucleo-synthesis, stellar evolution, neutrino-mass 
determinations and other fundamental processes of the Standard Model. 
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A Deriving the Faddeev Equation 

The Faddeev integral equation in the kinematics defined by Fig. ^2 is 

{t^''^Araa{E;k,p} = {K^^,A)'\a{E-Xv) (A.l) 

where E = — ^ is the total non-relativistic energy in the cm-frame, k (p) the momentum 
of the incoming (outgoing) deuteron/spin-zero field with spin-isospin indices i/A (j/B), and 
the incoming (outgoing) nucleon has spin-isospin aa {(3b). The integration over the loop- 
energy was already performed, setting Qo = 2m under the assumption that no additional 
poles above the real go-axis are hidden in the interaction K. Sub-scripts (super-scripts) 
denote quantum numbers for incoming (outgoing) particles. 

One decomposes now the integral equation into the pertinent spin-isospin and angular- 
momentum channels. 




Figure 11: The Faddeev equation for A^rf-scattering to N^LO. Thick solid line: propagator 
of the two intermediate auxiliary fields dg and d^, denoted by V, see fj3.6|) : triple line: 
propagator of the triton auxiliary field t, see ()3.17p : K: interaction to N^LO. 



A.l Projectors 

As two cluster-configurations exist, namely Ndt and Ndg, it is convenient to decompose 
each operator (interactions, amplitudes etc.) as 




(A.2) 
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In other words, each operator is represented in cluster-configuration space by a 2x2-matrix 
which in addition carries spin- and iso-spin indices. It is understood in the following that 
the spin-isospin indices on the 2x2-matrix are applied to each entry separately, and that 
Kronecker-(5's in spin or iso-spin are not displayed. Thus, for example, an entry 0{Ndt 
Ndt) = (JiCF^ in the above matrix is written out as 

{0{Ndt ^ NdtY^'^aa = {cy^Cy'fa 5^ (A.3) 

O is usually not symmetric in any pair of spin-isospin indices (zj), {AB), (a/5), {ah). 

The projection onto a state with angular momentum /, connecting momenta g,p with 
q-p = pq cos 6 and p = \p\ is as usual given by 

1 

0^^\q,p) = 1 J dcos^ 0{q,p) . (A.4) 
-1 

Finally, the projectors onto the possible spin-isospin states of the three-nucleon system 
are: combining the auxiliary fields with spin-index i or iso-spin index A with a nucleon of 
spin-isospin a, a into the spin-doublet channel with spin /x and iso-spin m 

1 /n- n \ "^^ 

(P...r^„^^(o J . (A.5) 
and into the spin-quartet channel with spin (fij) and iso-spin m 

The projectors are related to their Hermitean conjugates by 

= {V,m)^ = "P,,^" , {Vi,y = VI, , (A.7) 

and are ortho-normalised: 

V,,aV',^ = 1 , Vl^Vl = , V,,aVI, = (A.8) 

For a complete set of spin-isospin projectors, supplement these by the (ortho-normalised) 
projector onto the iso-spin quartet channel with iso-spin [mB) and spin /x, which is not 
found in nucleon-deuteron scattering: 

\ ^ 3 ^/ aa 
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A. 2 Projecting the Interaction Terms 

The exchange of a nucleoli is from the Lagrangeans ()2.1|/I2.2^2.2|) up to N^LO: 



N-exiA 



h0 



bp 



+ +p- q- ME -\e \a,T^ tat^ 



(A.IO) 



The projection of the propagator of the exchanged nucleon onto angular momentum / is 
obtained by combining ()A.4|) with (j3.2j) . There is no mixture between the doublet- and 
quartet-channels, and one finds in the doublet-channel 



V, 



(Ti(T^ Ta 

d,jB I ^ ^ 



-1 3 
3 -1 



(A.ll) 



and in the quartet 



^ a, tat'' ' ' 



(A.12) 



In cluster-configuration space, one obtains for the three-body interaction from ()3.16|) 



iB \ 

Kn iA ) {E]q,p) 



3Mn 



E 

m,n=l 



/^2n(A) 



{ME + Yt 



2\n 



a-' a.; 



-Oi T^ 



-crJ Ta 
r^TA 



bl3 



(A.13) 



By construction, the only non-zero contribution is in the doublet-S wave: 



V, 



-a^ Ta 



d,jB 



-a, T^ T^TA 



1 -1 
-1 1 



(A.14) 



A. 3 Result 

Putting these results together and multiplying in the doublet-case the 2 x 2-matrix t^'^ in 
cluster-configuration space from the right with the column vector , one projects finally 
onto the nucleon-deuteron system. There is no mixture or breaking between in- or out- 
states of different individual nucleon spin and iso-spin. One arrives thus finally at the 
integral equations for the quartet ()3.1|) and doublet ()3.4p channels quoted in the main text. 
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